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Abstract — We prove that if the non-self-adjoint scalar wave equation satisfies Huy- 
gens' principle on Petrov type III space-times, then it is equivalent to the conformally 



in 
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Or 



invariant scalar wave equation. 



Resume — On demontre que la valite du principe de Huygens pour l'equation des onde 
scalaires non-auto-adjointe sur un espace-temps general de type III de Petrov implique 
f*^ ■ que l'equation est equivalente a l'equation invariante conforme des ondes scalaires. 



1 Introduction 

In this paper we consider the general linear second-order hyperbolic equation on a 
■ four-dimensional curved space-time V4 , with C°° coefficients: 

Pu := g ab V a V b u + A a (x)d a u + C(x)u = , (1) 

where g ab is a contravariant pseudo-Riemannian metric with signature 
(-1 ) on V4, u is the unknown scalar function, V a denotes the covariant deriva- 
tive with respect to the Levi-Civita connection, A a is a vector field and C is a scalar 
field in V4. All considerations are restricted to a geodesically convex domain. The 
equation (JTJ) is also called the non-self-adjoint scalar wave equation. 

Huygens ' principle is said to be valid for an equation of the form (0) if and only 
if for every Cauchy initial problem, and for each point Xq G V4, the solution depends 
only on the Cauchy data in an arbitrarily small neighbourhood of S R C~(xo), where 
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S denotes the initial surface and C~(x ) the past null conoid of xq. Such equations 
are called Huy gens' equations. 

Necessary conditions for the validity of Huygens' principle for are given by 

(/) C :=C- l -A i . i - l l A i A i - l -R = Q, (2) 

(//) H k a , k = 0, (3) 

(III) S abk . k — -C k ab l L kl = —5 ^H ak H b k — -g ab H k iH klSs j , (4) 

(IV) TS(3S abk H k c + C k ab l H ck]l ) = , (5) 

(V) TS [?>C k ab l ™Ckcdi;m + 8C k ab l -cSkid + ^0S a b k S cdk — 8C k ab l S k i c -d 

riA/~yk In i a l/~im j i inril: I /~im j 

— ab Jcdk;l + ab L>1 ck^dm + ab (-> cdl-^km 

+ 12H ka . bc H k d — 16H ka . b H k c . d — 8AH k a C kbc iH l d 

-18H ka H k b L cd ) = 0, (6) 

(VI) TS (3QC k ab l Ci cdm]k H m e — 6C k ab l -cOide" 1 H km — 138S ab k C kcd iH l e 
+6S abk H k c . de + 6C k ab l - c Hkd;ie ~ ^S abk . c H k d . e 

+12C k a b L kc Hi d;e — 9C k a b l ;c L kd Hi e — 9S a b k L cd H k e ) = , (7) 



where 



Cijki '■— Rijki — 2g[i[iLj] k ] , (9) 

Sijk '■= Li[j- k ] , (10) 
R 

Lij := —Rij + — Qij . (11) 

In the expressions above R a bcd denotes the Riemann tensor, C abcd the Weyl tensor, 
Rab '■= g cd R ca bd, the Ricci tensor, and R := g ab R a b the Ricci scalar associated to the 
metric g ab , and A a := g a bA b . Conditions I — IV were obtained by Gunther [Zj, condi- 
tion V was derived by Wiinsch ^Sj for the self-adjoint case and McLenaghan |Hj for 
the non-self-adjoint case. Condition VI was obtained by Anderson and McLenaghan 
[2]. We shall use here the two-component spinor formalism of Penrose [Tl] and the 
spin-coefficient formalism of Newman and Penrose |13| ITo] . whose conventions we 
follow. The spinor equivalents of the tensors (jSJ), (JHJ), (jTUjl (jTTj) are given by : 

Cabcd «-> ^ABCD^AB^DC + ^ ABCD e AB^DC > (12) 

if a6 ^ 2((f} AB e Ai3 + \j> AB ) , (13) 
L ab <-> 2($ abAb - Ae AB e AB ) , (14) 

Sabc ^ D ABC;DA £ CB + * ABCjDA^CB , (15) 
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where ^abcd = ^(abcd) is the Weyl spinor, A = (1/24)R, (p AB = 4>{ab) is the 
Maxwell spinor, and &abAb = &(ab)(Ab) = ^abAb * s the trace-free Ricci spinor. 

We can set up, at each point of space-time, a dyad basis {o A , i A } satisfying 
the relation oa^ a = 1- The NP-components of any spinor are defined by projecting 
the spinor into the the local basis {o A , i A }. For the curvature spinors and Maxwell 
spinor we have 

^ABCD = ^0 L ABCD ~ ^lO(A^BCD) + ^20(AB^CD) 

-4:<S> 3 {A BCtD) + ^iOABCD , (16) 



®ABAB = ^220ABOj iB ~ 2$ ' 2 i0 ABO '^t B) ~ 2^ 12 O^B)0 AB 
+ ®20OabTab + ^02l'ABTAB + ^llO(A^B)0( A i B) 
-2$ 10 O(AtB)O( A l B) ~ 2§ 01 l A B'tAB > ( 17 ) 



Ub 



- 20!O(AiB) + 4>20AB , (18) 



where l A bcd = la^bLc^d, etc. 

The covariant derivatives of the dyad basis spinors are given in terms of the 
NP-spin coefficients by 



°A:BB - °a! BB + 1>aII B B > L A:BB ~ C 'A 111 BB ~ l A^ BB , (19) 



where 



I BB '■= lo B o B - ao B I B - (3l b o b + ei B T B , 

11 BB '■= - T °BO B + pO B I B + Ol B B - KL B T B , \ (20) 

III BB := vo B o B - \o B T B - [Li B o B + ra B I B . 

The necessary conditions I - VI can be expressed in terms of dyad components, 
containing only the Newman-Penrose scalars. This conversion procedure consists in 
two steps. Firstly, we have to convert the tensorial expressions into spinor form, 
using (fT2"j) - (fT5J) . Secondly, the spinor equations must be expressed in terms of dyad 
components, using (fTT)|) . (fTTj) IfTBjl . and We perform these lengthy calculations 

automatically with the NPspinor package jU E], available in the Maple computer 
algebra system. 

Using the fact that Huj.u = 0, it can be shown that the spinor form of condition 
// is given by: 

(lis) ( f )AK . K A = 0. (21) 
For condition /// a direct application of the correspondence relations yields 

^ABKL; A B + ^ABKL K A L B + ^AB KL ^ K LAB 

+^AB kL ^KLAB + 10<Wi* = . (22) 
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Instead of (j22|) we shall use a stronger form of this condition, obtained by Wiinsch 
[T9"] and McLenaghan and Williams [TTj : 

{Ills) V K A V L ^ABLK + ® KL AB^ABKL + ^Ab4>AB = ° ■ ( 2 3) 

While the original necessary condition f!22D is Hermitian, (|2Hjl is complex. 

The conversion of the remaining conditions to the respective spinor form and 
the determination of the dyad components is done automatically by defining templates 
in the NPspinor package. The spinor form of the trace-free symmetric part of a 
tensor is obtained by taking the correspondent spinor equivalent of that tensor and 
symmetrizing with respect to all dotted and undotted indices [T7j. McLenaghan and 
Walton [TOj have shown that any non-self-adjoint equation (JTJ) on any Petrov type 
N space-time satisfies Huygens' principle if and only if it is equivalent to a scalar 
equation with A a = and C = 0, on an space-time corresponding to the exact 
plane-wave metric. 

In this paper we prove the correspondent result for Petrov type III space-times. 
In this case the Weyl spinor has the form: 

^abcd = a( A aBacPD) , (24) 

where a a and (3a are the principal spinors. If we choose the spin basis such that a a 
is proportional to the dyad basis spinor oa and (3 proportional to l a , we obtain from 

^ABCD = -4&30(ABCLD) ■ (25) 

The dyad transformation 

o' = e w/2 o, t' = e - w/2 (L + qo), (26) 
where w and q are complex functions, induces the following transformation on \& 3 : 

% = e- w (^ 3 + 3q^ 2 + Sq 2 ^ 1 + q 3 ^o)- (27) 
Thus, we can choose the tetrad such that \&3 = — 1, so that 

^ABCD = ^O(ABC^D) ■ (28) 

In a previous paper, Anderson, McLenaghan and Walton [3] have proved the 
following theorems: 

Theorem 1 The validity of Huygens' principle for any non-self-adjoint scalar wave 
equation |7J) in any Petrov type III space-time implies that the space-time is confor- 
mally related to one in which every repeated null vector field of the Weyl tensor l a , is 
recurrent, i.e., 

l[Jb;c] = . (29) 
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Theorem 2 There exist no non- self- adjoint Huygens' equations (QJ) on any Petrov 
type III space-time for which the following conditions hold 

*ABCB;EE^ B l C D t E E = , ] 

V ABCD . iEE i A L B L C o D o E o E = 0, (30) 

= • J 

In the next section we show that the restrictions imposed by these two theorems can 
be removed. 

2 Main Theorem 

The main result of this paper is expressed by the theorem: 

Theorem 3 (Main Theorem) If a non- self- adjoint scalar wave equation of the 
form (CP satisfies Huygens ' principle on any Petrov type III space-time, then it must 
be equivalent to a conformally invariant scalar wave equation 

g ah V a V b u+\Ru = Q. (31) 
o 

Proof 

We shall first prove, using the necessary conditions II to VI given by Q - 
(J7J), that the assumption A^ := H ab ^ leads to a contradiction. In terms of the 
dyad components of the Maxwell tensor H ab , this is the same as proving that the 
necessary conditions imply <p Q = <pi = <p 2 = 0. Finally we invoke a lemma by Giinther 
[7j that states that every equation of the form (QJ) for which A^ b ] := is related by 
a trivial transformation to one for which A a = 0. It then follows from the necessary 
condition / (eq. ©) that C = R/6. We use a notation for the dyad components of the 
necessary conditions in the form X ab , where X is the Roman numeral corresponding 
to the necessary condition, a denotes the number of indices corresponding to the dyad 
spinor i and b the number of dotted indices corresponding to the dyad spinor I. We 
shall refer to the Newman- Penrose field equations using the notation NP1, NP2, etc. 
as listed in the Appendix. The methods employed in this proof are similar to those 
used in [12]. 

We start with the result obtained by Anderson and McLenaghan |Sj, expressed 
in following lemma : 

Lemma 1 For the non-self- adjoint scalar equation of the form (flj) ; in Petrov type 
III space-times, the necessary conditions II , III , IV , V and VI together with the 
assumption that the Maxwell spinor (pAB is nonzero, imply that there exists a spinor 
dyad {oa, i>a} and a conformal transformation such that 

K = a = p = T = e = 0, 

= #i = ^4 = 0, #3 = -1, 
$oo = $01 = $02 = A = , 
Da = Dp = m> n = . 
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In what follows we shall use the relations (}3*2*j) where necessary. It was shown in |3] 
that 0o = 0i = 0. Thus, what remains to be proved is that the assumption 2 7^ 
leads to a contradiction. Let us assume initially that a(37r 7^ 0. From IIqx, IIqo, IH10, 
(NP6) and (NP25) we have, respectively, 

D0 2 = , (33) 

<50 2 = -20 2 /?, (34) 

5(3 = -(3{a + (3) , (35) 

D 7 = am + (3tc + $ u , (36) 

8®2i = 2(a + 2ti + A$n - a$ 2 i) • (37) 

By adding (NP22) to the complex conjugate of (NP23), and solving for D$ 12 we get 

D$ 12 = 27f$ u , (38) 

and 

5$n=0. (39) 
Subtracting (NP24) from (NP29) and solving for 5$i 2 we obtain 

^i2 = 2(- / 9 + $n/Z-)9$i 2 ), (40) 
Adding (NP24) to two times (NP29) and solving for D$ 22 we get 

D$ 22 = 2(-/3-^ + $ 21 7f+$ 12 7r). (41) 
By substituting (jBH) into IViq we find 

50 2 = 2(3^0 2 - «0 2 + 2VT0 2 + «0 2 ) . (42) 

Using (f4*T]) and (fl^Jl . V 2 o and W03 can be written respectively as 

-240 2 2( 5 2 + 2 2 (-6mf + 18/3tt + 9 (3a - ea 2 - 2<5vf - 5a) 

+ 12a 2 + 2AWa + 80w(3 + +45a + 85w + U(3a = , (43) 

P{fa(8a + 25w + 3a 2 - 18W(3 + 6aW - 9a/?) + 240 2 /? 2 ) = . (44) 
Eliminating 5 (a + 2tt) from |43|) and (}44|) we find 

/3 2 — 

2 _ 4 (l9vr0 2 + lOa0 2 - 12/30 2 ) = . (45) 

The denominator — 2 2 + 4 in the expression above must be nonzero, since 2 cannot 
be constant. Otherwise, from (|34p. we would have (3 = 0. 
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In order to determine further side relations we still need to find the Pfaffians 
5a, 60, and 5n, in terms of 6a. From (NP6), (NP7), (NP8), (NP9) and (NP12) we 
have 

~5n = DA — 7r 2 — na + n~/3 , (46) 
6n = Dfi — nn + na — j3ir , (47) 
Du = An + n^i + n\i + 7f A + 7T7 - 7T7 - 1 + $21 , (48) 
Sa = 5P + aa + p]3- 2pa + $ n . (49) 

Using (J3~%|) and (jHJ), we now evaluate the NP commutator [<5, D]$ 22 — [A,D]$ 12 to 
obtain 

-25/3 + DA$i 2 + 27f$n7 + 20 + Ap 2 + 2ttA$h - QP$ 21 n + Ann 

+ 4/3a - D5$22 - 4o7$i 2 vr - 2/?$i 2 vr - 2a^ 2 {n + 27f$n7 - 27f$n/J 

+ 47f$n/x + 2¥$i 2 7r + 2$ 21 5¥ + 2$ 12 <57T - 2$n A¥ + 2a~]3 - 6n~/3 

+ 2$ 21 ¥ 2 + 27T$ 12 /3 + 2na = . (50) 

From [A, D]$i 2 we get 

DA$i2 = 2$n A¥ - 2¥$n7 - 2¥$ n 7 + 2na + Ann + 27rA$n 

-2m 12 n - 2nf3 + 2n$ ll Jl - 2n<S> l2 p . (51) 

From (NP26), 

5$ 22 = A$ 12 + 27 + 4/1 - 277$ n + 2A$ 2 i + 2$i 2 7 + 2$i 2/ u - 2$ 22 /3 - 2$ 22 a . (52) 

Substituting §2), flEEJ), (jHf|, (|47jl and (jiSjl into JSH we have 

Sp = -ap - An~P - 2DJI - (3p + 2nn - 2$ n . (53) 

Now, using (J33J), (jS3), (02), (EZ|), (03) and (J53J) in the commutator [6, 5}4> 2 , and solving 
for D/i we obtain 

D/i = ^(80 2 avf + 240 2( 5vr - 80 2 avf + lO0 2 $n + \2p~$ 2 a - 12a0 2 /3 

120 2 

+40 2 7T7r + 4o70 2 a - 4Q~</> 2 a - 2<fr 2 <&ii - 2A(fi 2 n~j3 + Scj) 2 na) . (54) 
The explicit form of DA can be determined from (|43|) . (|44j) and (|46|) . and is given by: 

DA = - J- fe 2 (3a 2 - 2tt 2 + Aan - I6f3n - 9pa + So) 

2(-0 2 2 + 4) v 

+8tt 2 - 16™ - 12a 2 + 240 2 2 /f - 88np - AApa - Ada)) . (55) 
Finally, from V\ 2 we get 

Dtt = 0. (56) 



We have now determined all the Pfaffians needed for finding new side relations 
using integrability conditions. Before proceeding we go back to ()45|) and introduce a 
further simplification by expressing 02 in terms of <j) 2 , 

- 2 =(19tt + 10a) . (57) 



12/3 

Substituting (JSTj) into the numerator of the complex conjugate of we find 

Si := 361vrvr + 190avr + 190vrct + lOOaa - 144/3/3 = . (58) 

Another side relation can be determined from the NP commutator [5 , 8](a + 2vr): 

(-820 a 2 $ n + 95vr 2 $n + 2680 a 3 a - 5448 a/3$ u - 10236 vr/3$n 
+ 5360 a 3 7f - 9648 P a 2 (3 - 38016 vr 2 /3 /3 - 13926 vr 2 /3 a - 2508 vr 2 vf /3 

- 3816 P a 2 a - 720 a 2 vr /3 + 61628 a vr 2 vr + 40128 vr 3 vr + 20064 vr 3 a 
+ 30814 vr 2 a a + 15752 tt a 2 a - 38376 vr P a p + 31504 tt a 2 vr 

- 14592 0^0"- 2688 vr a vf/3- 1508 vr a $ n ) /(19vr + 10a) =0. (59) 

It follows from (|57|) that the numerator 19 vr + 10 a in the preceding equation must 
be non-zero. Solving this equation for $n we obtain 

$n = (2680 a 3 a + 5360 a 3 vr- 9648/3 a 2 p- 38016 vr 2 /3/3- 13926 tt 2 /? a 

- 2508 tt 2 vf P - 3816 /3 a 2 a - 720 a 2 vf /3 + 61628 a vr 2 vf + 40128 tt 3 vf 
+ 20064 vr 3 a + 30814 vr 2 a a + 15752 tt a 2 a - 38376 tt /3 a p 

+ 31504 vra 2 vf- 14592avr/3a - 2688 vr a vf /3) /(820 a 2 - 95 vr 2 

+5448 /3a + 10236 vr/3 + 1508 vra), (60) 

where, for now, we assume that the denominator in the expression above, 

di := 820a 2 - 95vr 2 + 5448/fa + 10236vr/3 - 1508vra , (61) 

is non-zero. Evaluation of <50 2 + 20 2 /3 = (cf. Eq. PS])), using (jBTfl . and solving for 
$n gives 

$ n = (700a 2 a + 5300avrvf + 2650vraa + 2508vr 2 a + 1400a 2 vf 

- 2520 Pa p - 4752 vr /3 /3 - 1650 vr /3 a - 132 vr vf/3 + 5016 tt 2 vr 

- 900/3aa-72avf/3) /( -437 vr + 372/3 - 230 a ) , (62) 

where we assume, for now, that the denominator of ()62|). given by 

d 2 ■= 437vr - 372/3 + 230a ^ , (63) 
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is nonzero. Evaluating the commutators [5 , 8\(3 and [5 , 8\<f> 2 , and solving each one 
for 5 a we get, respectively, 

5a = (-308 7T a $n + 8016vr a 2 ¥+ 234 /3a 2 /3 + 3972 vr a 2 a + 78 /3a 2 a 
+ 759a/?$n - 29a 2 $n + 1440a 3 vr + 702a 3 a + 1386vr/3$n 
- 513 tt 2 $n + 1296 vr p a (3 + 528 vr 2 p a + 14872 a vr 2 If + 7436 vr 2 a a 
+ 1584vr 2 /3;# + 4598vr 3 a + 9196vr 3 r + 432avr/3a - 108vra 2 /3 
-54 a 3 12 vr a + 6 a a + 36 /? vr + 18 /3 a + 3 $ u ) , (64) 

— 6a 2 vr- 572vrvr^- 2867r^a-157/3aa + 3a 2 a-314avr^ 

a + 2 vr 

where we assume, for the moment, that the denominators of (J5l"j) and ([65)1 . given by 

d 3 := 4vra + 2aa + 12/3vr + 6/3a + $n , (66) 
<i 4 := a + 2vr , (67) 

are nonzero. By subtracting ()64|) from (|65|) and solving for $n we have 

$ n = (-864 /3a 2 a + 7436vr 2 aa - 3168avr/3a + 4008vra 2 a 

- 2904 vr 2 ]3 a + 720 a 3 a + 9196 vr 3 vr - 8712 vr 2 (3 ]3 + 14872 a vr 2 vr 
+ 4598 vr 3 a - 2592 P a 2 (3 + 80 16 vr a 2 vr + 1440 a 3 vr - 9504 vr ]3 a 0) 
/ /(-288/?a + 20a 2 - 528vr/? + 513vr 2 + 308vra) , (68) 

where the denominator of 

d 5 := -288/?a + 20a 2 - 528vr/3 + 513vr 2 + 308vra, (69) 

is assumed to be nonzero for the moment. Subtracting ()60|) from ()62|) and taking the 
numerator we find 

S 2 : = 39914208 vr 2 a 2 vf + 19957104 vr 2 a 2 a + 35332704 vr 3 vf/3 

+ 12279168 vr 3 /3 a + 1190400 a a 4 - 12773376 ~f a vr a - 3483648 p 2 a a 2 
+ 1200960 /?aa 3 - 11708928 /f a vr 2 - 20739456 vr/? a 2 (3 
+ 8008704 a 2 vr /? a + 30335616 vr a 2 vr J3 + 56708640 vr 2 a Wp 

- 32440608 vr 2 /3a p - 16161552 vr 3 pp + 8022720 vr a 3 a + 5408640 a 3 vf p 
+ 16045440 vr a 3 vr + 8529708 vr 4 a + 43221504 a vr 3 vr + 21610752 vr 3 a a 

+ 17059416 vr 4 vr + 2380800 a 4 vr + 17343792 avr 2 /3a - 10139904 ~f a 2 p 

- 124416 a 2 vf/f -418176 vr 2 vf/f -456192 vr a vf/f -4285440 pa 3 p 

- 37407744 vr p 2 a p - 34499520 vr 2 p p 2 = . (70) 
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By subtracting (p)2"j) from ()68|) and taking the numerator, we obtain a third side 
relation: 

S 3 ■ = -9374472 vr 2 a 2 W - 4687236 vr 2 a 2 a + 6137076 vr 3 tt P + 5150178 vr 3 /? a 

- 179600 a a 4 - 2128896 /J 2 ana - 580608 /f a a 2 + 686160/3 a a 3 

- 1951488 a vr 2 + 4189824 np a 2 (3 + 4040184 a 2 np a 
+ 5272368 na 2 n]3 + 9852984 vr 2 a vf /3 + 8913384 vr 2 /? a /3 

+ 6244920 tt 3 p ]3 - 1505080 tt a 3 a + 940320 a 3 n p - 3010160 vr a 3 vf 

- 3295930 n A a - 12877972 a tt 3 tt - 6438986 vr 3 a a - 6591860 vr 4 n 

- 359200 a 4 n + 7909932 a vr 2 /? a - 1689984 p 2 a 2 P - 20736 a 2 vf /f 

- 69696 7t 2 Wp 2 - 76032 tt a vr ]f + 646560 p a 3 P - 6234624 vr /f a /3 

- 5749920 tt 2 /? /f = 0. (71) 

We can eliminate n in the above equations by defining new variables x x and x 2 by 

xi:=-, (72) 

7T 

x 2 :==. (73) 

7T 

The side relations now assume the form (modulo non-zero factors) 

S x := 361 + 190xi + 190x7+ 100xi XT - 144x 2 XT = 0, (74) 

S 2 ■= -1663092x^x7- 99200 xlx 4 - 2527968 x 2 xa"+ 975744 X2" 2 xT 

- 1023264 XTxT + 2874960 x 2 - 4725720 x x x^ - 198400 x\ 

- 3326184 x\ - 1337120 x\ + 34848 XT 2 - 3601792 x x - 710809 XT 

- 2944392 xT + 844992 xT 2 x 2 x 2 + 1064448 xT 2 xTxi 

+ 290304 x? xix\ - 100080 x^xTx 3 + 1728288 XT x 2 x 2 

- 667392 x\xi XT + 2703384 xTxi x 2 + 1346796 x 2 XT — 668560 x\ XT 

- 450720 x\ xT - 1800896 x x x7 + 10368 x\ xT 2 + 38016 xi xT 2 

- 1445316 xixTx7 + 357120 XT x^x 2 

+ 31 17312 xi 2 xix 2 - 1421618 = 0, (75) 

5 3 := 2343618 x 2 x7+ 89800 x7x 4 - 2636184 x\x^ + 975744 xT 2 x7 

- 2575089 XTxT + 2874960 x 2 XT 2 + 3295930 - 4926492 x x XT + 179600 x\ 
+ 4687236 x\ + 1505080 x? + 34848 XT 2 + 6438986 x x + 1647965x7 

- 3068538 xT + 844992 xf x\ x 2 + 1064448 xf xl x x 

+ 290304 XT 2 xlx\ - 343080 x^xlx\ - 2094912 XT x\x 2 
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- 2020092 x\xl xi- 4456692 xixix 2 - 3122460 x 2 xi + 752540 x\xl 

- 470160 x\x^ + 3219493 x 1 xj + 10368 x\ xf + 38016 x x xf 

- 3954966 x x x^xj - 323280 x^x\ x 2 + 3117312 xf x x x 2 = . (76) 

In order to study the possible solutions of the polynomial systems that appear from 
the side relations, we apply the procedure gsolve which is part of the Maple package 
grobner j^J. This procedure computes a collection of reduced lexicographic Grobner 
bases corresponding to a set of polynomials. The system corresponding to the set 
is first subdivided by factorization. Then a variant of Buchberger's algorithm which 
factors all intermediate results is applied to each subsystem. The result is a list of 
reduced subsystems whose roots are those of the original system, but whose vari- 
ables have been successively eliminated and separated as far as possible. This means 
that instead of trying to find a Grobner basis, the package attempts to factor the 
polynomials that form the system after each step of the reduction algorithm. In the 
algorithm, the variables x\, x 2 and their complex conjugates, x~i and xi, are treated 
as independent variables. In the subsequent analysis we use the fact that they are 
complex conjugates of each other. 

Applying gsolve to the set of equations formed by S 2 = 0, S3 = (cf. (J75j) 
and ()76|)). their complex conjugates, and Si — (cf. Q74|)). we find the the only 
possible solution for which x\ 7^ and x 2 7^ is given by 

324x2^-1 = 0, 6xi + ll = 0, 6x1+11 = 0. (77) 

By substituting (|77jl into any of the previous expressions for $n we find that $n = 0. 
Using this and 7r = —a 6/11 in Vn one gets 

- 1761/3aa + 5a 2 a + 3267/3/f - 5445/3a/? = . (78) 

It is easy to verify that ((75)) and the first equation in (|77|). now given in the form 
1089/3/3 -aa = 0, imply that a = = 0. 

Let us consider first the cases in which each one of the denominators di, d 2 , 
g?3, di and ofe, given respectively by (fBTjl . (jHEJ), (jHTJ) and (jEHJ), is zero. 

(i) d x = 

From (j6*T|) we have, in terms of the variables x x and x 2 . 

820 x\ + 1508 X! + 5448 x^xi - 95 + 10236 x^ = . (79) 
Since the numerator of (|60|) must also vanish, we have 

2680 xf + 1340 x\ xj- 1908 x^xlxl + 7876 x\xl 

- 360x2^ + 15752 x\ - 4824x2"x 2 x\ + 30814 x x - 1344 x^X! 

+ 15407xi xT- 19188x^x 2 xi - 7296 xaxTxi + 20064 + 10032x1 

- 19008xix 2 - 6963 xi xi- 1254 xi = 0. (80) 
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Applying gsolve to the set of equations consisting of (|79j) . (JHUj) . their complex conju- 
gates, and S*i = (cf. eq. (JZH)) we find that all possible solutions require that x 2 = 0. 

(ii) d 2 = 

In this case, from (|Fk3|) ). we have 

230xi + 437-372x^=0. (81) 

This implies that the numerator of (|fi2*|) must vanish, so we have 

350 x\ xi + 700 x\ + 1325 x x xl - 36 xix x + 2650 x x + 2508 
— 45OX2X1X1 — 1260x"2X 2 Xi — 825x2x7 — 66x2 — 2376x2X 2 
+ 1254X1=0. (82) 

Applying gsolve to the set of equations consisting of (}8~Tj) . (J%2~j) . their complex con- 
jugates, and Si = we find again that all solutions require that x 2 = 0. 

(iii) d 3 = 

Eq. (JBT)|) now gives 

$u = -27T7r(xi + 2) (3x 2 + xl) . (83) 

By subtracting (jH3^1 from its complex conjugate we obtain 

Ei := —6x2 — 3xiX2 — 2xT + 6x2 + 3xix7 + 2xi = . (84) 

Subtracting (J%3~j) from (j6TI|) and taking the numerator, we get 

E 2 := 246x 2 Xi + 268xi + 216xiX7 2 + 477xiX 2 + 1152x2"XiX 2 
—36x2X1 + 1066xi + 354x2Xix7 + 2232x2X2 + 518x7 — 66x2 
+711x^x1- 30x 2 + 692xix7+ 1056 = 0. (85) 

An additional side relation is obtained by subtracting the complex conjugate of 
from (p)2*j) and taking the numerator of the resulting expression: 

E 3 := 144x 2 X2" + 2691xix 2 + 2622x 2 + 144x2XiX 2 + 690x 2 x? 
-120x^x7 - 428xix7 + 81x 2 Xi - 380x7 - 2508 + 66x^ - 2650xi 
+36x2-xi - 700x1 + 78x 2 xix7 = . (86) 

Applying grobner to the set consisting of E 2 , its complex conjugate, E x , E 3 and 5*1 
we find that this system admits no solution. 

(iv) d 4 = 



12 



When the denominator of ()65|). given by (joT|) . is zero, its numerator must be 
zero, implying that a = /3(443/3). This, on the other hand implies immediately, from 
Si = (cf. (|HH|), that p = 0. 

(v) d 5 = 

When the denominator of (|68|) is zero, its numerator must be zero too. Thus, 

we get 

(11 + 6xi) 2 (38 + 19xT'12xixT + IOxixT + 20xi - 36x 2 xi) = , (87) 

308xi + 513 + 20x1 - 528x^ - 288x^Xi = . (88) 

Applying gsolve to the polynomial system defined by the system of polynomials 
defined by (|57j). (fSSj). their complex conjugates, and Si = 0, we find that there are 
no possible solutions. 

The case a = (3 = tt = 0, which leads to Theorem |21 was considered in |3] and 
results in a contradiction to the hypothesis (f>2 7^ . The more general case, air(3 = 0, 
also leads to a contradiction. The proof, found in p]Gj|, is tedious but straightforward, 
and will not be presented here. Thus, the necessary conditions / — VI for the validity 
of Huygens' principle imply that we must have H a b := = 0. 

The last step of the proof requires the use of the following lemma [7J : 

Lemma 2 Every scalar wave equation of the form (QJ) for which A^j] = 0, is related 
by a trivial transformation to one for which Ai — 0. 

From necessary condition I (cf. eq. (J2J) we now have B = R/6, so the wave 
equation is conformally invariant. Thus, the Main Theorem is proved. 
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Appendix 

In this Appendix we give the Newman-Penrose field equations and commu- 
tation relations referred along the paper. We note that many lists in the literature 
contain typographic mistakes, or use different conventions. 

Bianchi identities 

(NP1) Dp - 5k = p 2 + aa + (e + e)p - kt - (3a + /3 - tt)k + $ 00 , 
(NP2) Da~5n= (p + p)<j + (3e - t)a - (r - W + a + 30)k + * , 
(NP3) Dr - Ak = (t + W)p + (r + n)a + (e - e)r - (37 + t")k + ^1 + $01, 
(AP4) Da - 5e = (p + e - 2e)a + /5a - /3e - kA - K7 + (e + p)it + $10, 
(NPb) Dp - 5e = (a + 7r)<7 + (p - e)(3 - (p + 7)^ + (tt - a)e + #1, 
(AP6) D7 - Ae = (r + 7f)a + (r + 7r)/3 - (e + e)7 - (7 + T)e + rvr - i/k 
+^ 2 _- A + $11, 

(NP7) DX- 5ti = pX + Wp + n 2 + (a - fi)ir -vK+(e- 3e)A + $ 2 o, 
(NP8) Bp - 5n = pp + aX + vrvf - (e + e)p - (a - /3)tt - vk + ^2 + 2A, 
(AP9) Dis - An = (r + 7r)/i + (r + w)X + (7 - j)-n - (3e + e)i/ + ^3 + $21, 
(AP10) AA -_<Si/ = (7 - 37 - p - p)A + (3a + (3 + tt - r)z/ - tf 4 , 
(AP11) 5p - 5tr = (a + /3)p - (3a - ~fi)o + (p - p)r + (p - p)n - *i + $01, 
(AP12) 5a-5(3 = pp- aX + aa + (3]3 - 2a[3 + (p - ph + (/i - p)e 
-^2 + A + $n, 

(AP13) 5X-6fjL = (p- J))u_+ (p - p)tt + (a + (3)p + (a - 3/3) A - ^3 + $21, 
(AP14) Su - Ap = p 2 + XX + (7 + 7)/! - Vn + (r - a - 3/3)// + $ 22 , 
(AP15) ^7 - A/3 = (r - a - (3)-f + pr - av - eV - (7 - 7 - p)/3 + a A + $ i2 , 
(AP16) 5r - Ao- = per + Ap + (r - a + /3)r - (37 - 7)0- - kT> + $ 02 , 
(AP17) Ap - 5t = -pp - aX + (7 + 7)p - (r + a - /3)t + i/k - ^ 2 - 2A, 
(AP18) Aa - 5 7 = (e + p)v - (r + /3)A + (7 - p)a + (/3 - r) 7 - ^ 3 • 



Ricci identities 

(AP19) 5^ - D^i + D$ i - 5$oo = (4a - tt)^ - 2(2p + e)*i 

+3k^ 2 + (w - 2a - 2/3)$ 00 + 2(e + p)$ i + 2a$i - 2/t$ n - k% 2 , 

(AP20) A^ - + D$ 02 - 5$oi = (47 - p)^ - 2(2r + 

3(7^2 - A$ 00 + 2(7f - /3)$ i + 2<7$ n + (2e - 2e + p)$ 02 - 2k$ 12 , 

(AP21) 35^i - 3D^ 2 + 2D$ n - 25<5> l0 + 5$oi - A$ 00 = 3A^ - 9p^ 2 

+6(a - tt)*i + 6^3 + (p - 2p - 2 7 - 27)$ 00 + (2a + 2tt + 2r)$ i 
+2(r - 2a + ¥)$ 10 + 2(2p - p)$ n + 2ct$ 20 - cf$ 2 - 2k$ 12 - 2k$ 2 i , 

(AP22) 3A^! - 35^2 + 2D$i2 - 25® u + M 02 - A$ i = 3z/^ + 6(7 - p)*i 
-9r^ 2 + 6cr^3 - ^$00 + 2(p - p - 7) $„i - 2A$ 10 + 2(r + 2w)$ u 
+ (2a + 2tt + t - 2/3) $02 + (2p - 2p - 4e)$i 2 + 2a$ 2 i - 2k<5> 22 , 
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(NP23) 3M 2 - 3D^ 3 + D$2i - <5$ 20 + 25$ n - 2A$ 10 = 6A*i - 9tt^ 2 

+6(e - p)* 3 + 3/t^4 - 2z/$ 00 + 2(p - p - 27)$i + (2tt + 4r)$ n 
+ (2/3 + 2r + vf - 2«)$ 20 - 2a$i2 + 2(p - p - e)$ 2 i - 7t$ 22 + 2A$ i, 

(NP24) 3A^ 2 - 35^3 + D$ 22 - 5§ 2 i + 25$i 2 - 2A$n = 6i/*i - 9p* 2 

+6(/3 - r)* 3 + 3cr^ 4 - 2z/$ 01 - 2z7$ 10 + 2(2/1 - p)$ n + 2A$ 02 - A$ 20 
+2(tt + r - 2/3)$i2 + 2((3 + r + 7f)$ 2 i + (p - 2e - 2e - 2p)$ 22 , 

(NP25) _ D* 4 + 5$2i - A$ 20 = 3A* 2 - 2(a + 2tt)* 3 + (4e - p)* 4 
-2i/$i + 2A$n + (27 - 27 + p)$ 20 + 2(T - a)$2i - 0^22 , 

(NP26) A* 3 - 5*4 + 5$ 22 - A$2i = 3z/* 2 - 2(7 + 2p)* 3 + (4/3 - r)* 4 
-2z/$n - I7$2o + 2A$i2 + 2(7 + p)$ 2 i + (r - 2/3 - 2a)$ 22 , 

(NP27) D$ n - 5$ 10 - 5$ i + A$ 00 + 3DA = (2 7 - p + 27 - p)$ 00 

+ (tt - 2a - 2r)$ i + (tt - 2a - 2r)$i + 2(p + p)$ n + a<l>02 + o-$ 20 

-7t$i2 - K$21 , 

(iVP28) D$ 12 - 5$n - 5$ 2 + A$ i + 35A =_(2 7 - p - 2p)$ i+ 
I7$ 00 - A$ 10 + 2(tt - r)$ n + (tt + 2/3 - 2a - r)$ 02 
+ (2p + p - 2e)$ 12 + 0"$2l - k$ 22 , 

(iVP29) D$ 22 - 5$ 21 - 5$i2 + A$ n + 3AA = i/$ 01 + 77$ 10 - 2(p + p) 
$ n - A$ 02 - A$ 20 + (2tt - T + 2^)$ 12 + (2(3 - T + 2tT)$21 
+ (p + p- 2e-2e)$ 22 . 



NP commutation relations 

55 -55 = (-p + p)D + (-p + p) A + (a - /3)<5 + (-5 + (3)5 , 
5 A - A5 = -uD + (t - a - /3) A + \5 + (p + 7 - 7)6 , 

- D5 = (a + /3 - tt)D + TtA - a5 - (p - e + e)5 , 
AD - DA = (7 + 7)D + (e + e) A - (r + n)5 - (r + 7f)£ . 
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